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ABSTRACT 

We have found a new clue to the transition mechanism between optical high/X-ray off and optical 
low/X-ray on states of the LMC supersoft X-ray source RX J0513.9— 6951. A sharp ~ 1 mag drop is 
common to the CI Aql 2000 outburst. These drops are naturally attributed to cessation of optically 
thick winds on white dwarfs. A detailed light-curve analysis of CI Aql indicates that the size of a disk 
drastically shrinks when the wind stops. This causes ~ 1 — 2 mag drop in the optical light curve. In RX 
J0513.9— 6951, the same mechanism reproduces sharp ~ 1 mag drop from optical high to low states. We 
predict this mechanism also works on the transition from low to high states. Interaction between the 
wind and the companion star attenuates the mass transfer and drives full cycles of low and high states. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: individual (CI Aquilae, 
RX J0513.9— 6951) — stars: winds, outflows — X-rays: stars 



1. INTRODUCTION 

The Large Magellanic Cloud (LMC) supersoft X-ray 
source RX J0513.9-6951 (hereafter RX J0513) shows a 
prominent recurrency of optical high and low states with 
quasi-regular intervals (~ 120 days high and ~ 40 days 
low states). Copious supersoft X-rays were detected only 
in the optical low states (e.g., Schaeidt et al. 1993). 
The transition mechanism between high and low states 
have not been fully elucidated yet, although a few ideas 
were proposed (see, e.g., Reinsch et al. 2000, for re- 
cent progress). We have found a new clue to this tran- 
sition mechanism from modeling of the recurrent nova CI 
Aquilae. 

CI Aql erupted in 2000 April has been densely observed 
in various optical bands (e.g., Hachisu, Kato, & Schaefer 
2002; Kiss et al. 2001; Lederle & Kimeswenger 2002; Mat- 
sumoto et al. 2001; Schmeja, Armsdorfer, & Kimeswenger 
2000) and in X-ray bands (Greiner k, DiStefano 2002). It 
reached the optical maximum (V ~ 9 mag) on 2000 May 
5. The visual brightness quickly decreased to V ~ 13.5 in 
about 50 days. Then, it stays at V ~ 14 for about 150 
days, i.e., a plateau phase. 

Matsumoto et al. (2001) reported a sharp ~ 2 mag drop 
of i? c -magnitude on 2000 November 23, about 200 days 
after the optical maximum, and Kiyota (2001, VSNET 
archives, http://vsnet.kusastro.kyoto-u.ac.jp/vsnet/) also 
observed a ~ 2 mag drop of J c -magnitude around the same 
day as Matsumoto et al.'s (2001) observation (see, e.g., 
Hachisu et al. 2002). Shortly after this drop, both the 
R c - and 7 c -magnitudes recovered by ~ 1 mag but stayed at 
~ 1 mag below the level before the drop. Hachisu & Kato 
(2001a) and Hachisu et al. (2002) attributed this drop to 



cessation of massive winds by calculating time-evolution 
of optically thick winds on a 1.2 M Q white dwarf (WD). 

The above feature of CI Aql light curve reminds us the 
optical high-to-low transition of RX J0513. In this paper, 
we try to reproduce the sharp drop in the optical light 
curves, both for CI Aql and RX J0513, by the same mech- 
anism. In §2, our numerical method and models for the CI 
Aql 2000 outburst are briefly introduced and summarized. 
The numerical results for the CI Aql 2000 outburst are 
given in §3. We try to reproduce the transition between 
high and low states of RX J0513 in §4. Discussion follows 
in §5. 

2. BASIC MODEL OF CI AQUILAE 2000 OUTBURST 

The mid plateau phase is a common feature among the 
U Sco subclass of the recurrent nova class, to which both 
U Sco and CI Aql belong. The plateau phase is well ex- 
plained by the contribution of an irradiated disk (e.g., 
Hachisu et al. 2000). The presence of the disk in out- 
burst was observationally confirmed by Thoroughgood et 
al. (2001). 

Our binary model is illustrated in Figure 1, which con- 
sists of a main-sequence star (MS) filling its Roche lobe, 
a WD photosphere, and a disk around the WD. A circu- 
lar orbit is assumed. We also assume that the surfaces of 
the WD, the companion, and the disk emit photons as a 
blackbody at a local temperature which varies with posi- 
tion. We assume an axi-symmetric accretion disk with the 
size and thickness of 

-Rdisk = aR*, (1) 

and 

h = /3i?disk ( -^-) , (2) 
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where i?disk is the outer edge of the optically thick part of 
the accretion disk, R\ is the effective radius of the inner 
critical Roche lobe for the WD component, h is the height 
of the surface from the equatorial plane, and w is the dis- 
tance from the symmetry axis. Here, we adopt v = 1 
during the strong wind phase, but v = 2 when rapid mass 
accretion resumes after the optically thick wind stops, to 
mimic the flaring- up disk by spray (e.g., Schandl, Meyer- 
Hofmeister, & Meyer 1997). The other two parameters of 
a and [3 arc determined by light curve fittings. 

We have used the same parameters as those in Hachisu 
et al. (2002) unless otherwise specified. The mass of 
the white dwarf was determined to be 1.2 ± 0.05 Mq by 
Hachisu & Kato (2001a). The mass and temperature 
of the companion star were estimated to be ~ 1.5M Q 
and 7,300 K (see, e.g., Hachisu & Kato 2001b, for the 
present evolutional status) from the light curve fitting of 
the orbital modulations in the 2000 outburst (Hachisu et 
al. 2002; Lcderle & Kimeswenger 2002). The inclination 
angle is a free parameter, which is redetermined from our 
light curve fitting. 

Time evolution of the photospheric radius and temper- 
ature of the WD are calculated from Kato & Hachisu's 
(1994) optically thick wind solutions (see Hachisu & Kato 
2001b, for more details of the numerical methods). When 
the photospheric radius shrinks smaller than the binary 
size, the irradiation effects both of the disk and of the 
companion become important (see, e.g., Fig. 2 of Hachisu 
et al. 2000, for such an example). 

3. WHAT WE LEARNED FROM CI AQUILAE 

Orbital modulations during the plateau phase have been 
reported by Matsumoto et al. (2001). Their R c light curve 
data show a scatter of ~ 0.1 — 0.2 mag so that we aver- 
age them into 36 phase bins as shown in Figure 2. No 
secondary eclipses were observed in contrast to the later 
phase orbital light curves (Lederle & Kimeswenger 2002). 
To fit these orbital modulations as well as the entire evolu- 
tion of i? c -magnitude, we have calculated light curves with 
various parameters (a is changed by 0.1 steps, (3 by 0.05 
steps, and the inclination angle i by 1° steps). The best 
fit model is shown in Figure 2, i.e., (a, (i) — (2.0, 0.05) and 
i = 69°. 

Here we plot two identical orbital modulations, but one 
is shifted by 0.36 mag down. The brighter one corresponds 
to the earlier stage on 2000 July 17 while the fainter one 
does to that on 2000 October 1. To reproduce the light 
curves with no secondary eclipses, we need both (1) large 
(area) irradiation effect to reproduce the brightness, and 
(2) bright sources other than the companion to eliminate 
the secondary eclipse. Then, we assumed that (a) the outer 
edge of an optically thick disk extends up to almost the 
orbital separation as illustrated in Figure la, and that (b) 
the WD photosphere is as large as i?wD, P h ~ 0.35 Rq 
on July 17 and i?wD,ph ~ 0.15 Rq on October 1 if the 
other conditions are the same between these two epochs. 
Any light sources other than the disk or the companion 
are unlikely to contribute to the i? c -magnitudc because 
the strong wind blows up off nebulosities around. This 
large extension of the disk naturally explains evident lack 
of secondary minima. The width of sharp eclipse minima, 
~ 0.05 — 0.1 orbital phase, is caused by the partial occul- 



tation of the WD photosphere. These large photospheric 
radii (i?wD, P h ~ 0.35 — ► 0.15 Rq from July 17 to October 
1) are very consistent with our theoretical model because 
WDs blow a wind when i?wD, P h > 0.085 Rq (e.g., Kato 
& Hachisu 1994; Hachisu & Kato 2001b; Hachisu et al. 
2002). 

The effect of a larger inclination angle is canceled out 
by the effect of a larger radius of the disk. For example, a 
combination of the inclination angle i = 74° and the disk 
of (ot,@) = (3.0,0.05) can also reproduce the light curve 
modulations (dashed line) as shown in Figure 2. Its eclipse 
dip seems to be too deep to be compatible with the obser- 
vation. However, the total amplitude of the modulations 
is ~ 0.8 mag (see Fig. 4 of Matsumoto et al. 2001, for the 
original observation) and still consistent. Here we adopted 
a set of i = 74° and (a, f3) = (3.0,0.05). In both i = 69° 
and i = 74° cases, the secondary occults only a part of the 
disk. This is, however, enough to make a orbital modula- 
tion observed because the central part of the disk is much 
brighter than the outer part. 

The velocities of winds calculated are as fast as 1000 — 
1500 km s _1 (e.g., Kato & Hachisu 1994) and the wind 
mass loss rate is as large as 10~ 6 — 1O~ 4 M yr~ x (e.g., 
Hachisu & Kato 2001b). Observationally much faster 
wind velocities are reported, for example, 2000 — 2500 
km s" 1 in CI Aql (see, e.g., Kiss et al. 2001). Such a 
strong wind affects the surface of the accretion disk. Be- 
cause of a large velocity difference between the wind and 
the disk surface, it certainly drives the Kclvin-Helmholtz 
instability at the interface. Only the very surface layer 
of the disk is dragged away outward with the velocity at 
least several hundred km s _1 like a free stream moving 
outward. This surface free stream is optically thick near 
the original disk but becomes optically thin somewhere 
outside because of geometrical dilution effect. We regard 
the transition place from optically thick to thin as the edge 
of the extended disk. 

It should be noted that a high density part of the 
disk is hardly changed by this relatively tenuous surface 
free stream and still resides within its Roche lobe be- 
cause the internal density of the disk is much denser than 
that of the WD wind. The wind mass loss rate is about 
1 x 10~ 6 M Q yr _1 and its velocity is ~ 1000 km s _1 , so 
that the density of the wind is estimated from the conti- 
nuity (M w i n( j = 4irr 2 pv) to be about 10~ n g cm 3 at the 
distance of 1 Rq from the center of the WD. On the other 
hand, the density of the standard accretion disk is about 
1 x 10 _3 g cm 3 at the same radius. Here, we assume the 
standard accretion disk model (Shakura & Sunyaev 1973) 
with a mass accretion rate of ~ 1 x 10~ 7 -M© yr _1 . 

After the wind stops, the photospheric radius of the WD 
quickly shrinks in the Kelvin-Helmholtz timescale of the 
envelope, i.e., 

A£ onv ,thormai 2.3 x 10 44 erg 

TKH = r ~ o T7TS T ~ 13 da y S ' ( 3 ) 

L 2 x 10 38 erg s^ 1 

where tkh and Ai? cnVjt hcrmai are the Kclvin-Helmholtz 
timescale and the thermal energy of the WD envelope, 
respectively, and L is the luminosity of the WD. Here, we 
adopt a typical thermal energy of the WD envelope from 
our solutions. It takes 20 days from i?wD,ph = 0.085 Rq 
to 0.04 Rq and takes another 20 days from 0.04 Rq to 
0.02 Rq from our wind solutions and evolution models. 
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The disk edge of optically thick part shrinks to a normal 
size of 0.7 — 0.8 times the Roche lobe radius as observed 
in the supersoft X-ray sources (e.g., Hachisu et al. 2002; 
Lederle & Kimeswenger 2002). This may occur in several 
dynamical timescales, i.e., a few to several days. Here, 
we adopt 10 days, i.e., the transition from configuration 
(a) to (b) in Figure 1. We adopt {a, (3) = (0.8,0.05) for 
configuration (b). Then, the mass transfer rate from the 
MS to the WD increases to make a spray around the ac- 
cretion disk as seen in luminous supersoft X-ray sources. 
This can be mimicked by a flaring edge of the accretion 
disk, i.e., configuration (c) in Figure 1 (see, e.g., Hachisu 
et al. 2002; Lederle & Kimeswenger 2002, for light 
curves and flaring- up disk models after the wind stops). 
We adopt (a, [3) ~ (0.8, 0.3). It also takes several dynami- 
cal timescales. Here we take 10 — 20 days for the formation 
time of the spray. Thus, we are able to reproduce a sharp 
~ 2 mag drop around 2000 November 23 and the ensuing 
~ 1 mag recovery soon after the drop, as shown in Figure 
3. A similar transition was observed in the U Sco 1999 
outburst on HJD 2,451,244 just before the detection of su- 
persoft X-rays (see Fig. 2 of Hachisu et al. 2000; Kahabka 
et al. 1999). 

4. LMC SUPERSOFT X-RAY SOURCE RX J0513. 9-6951 

The LMC supersoft X-ray source RX J0513. 9-6951 
shows regular transitions between 60 — 120 days optical 
high (X-ray off) and 30 — 40 days optical low (X-ray on) 
state. The X-ray off/optical high state was understood by 
a large expansion of the WD photosphere, which is trig- 
gered by an increase in the mass transfer rate onto the WD 
(e.g., Meyer-Hofmeister et al. 1997; Pakull et al. 1993; 
Reinsch et al. 1996, 2000; Southwell et al. 1996). When 
the WD envelope expands, an optically thick wind begins 
to blow. As we already know in CI Aql, the sharp drop 
of V-m&g is naturally explained by cessation of the wind. 
Our scenario is as follows: 

1. We assume the metallicity of Z = 0.004 and hydrogen 
content of X = 0.7 for the envelope of the WD in RX J0513 
since the LMC metallicity is reported to be about a third 
of the solar metallicity. 

2. We have directly calculated the duration of the wind 
phase and the Kelvin-Helmholtz timescale on massive 
WDs of 1.0, 1.1, 1.2, 1.3, 1.35, 1.36, 1.37, and 1.377 M 
(see also Hachisu & Kato 2001b, for wind solutions). The 
1.3 M Q WD model gives us a reasonable wind timescale 
of T w ; nc i ^120 days and a Kelvin-Helmholtz timescale of 
tkh ~ 13 days. 

3. We use the ephemeris revised by Cowley et al. (2002). 
The inclination angle is not known but suggested to be 
low. Here we adopt i = 10°. The light curve is not 
so different among i = 0° — 20°. The calculated light 
curve is plotted in Figure 4, which is made by connect- 
ing brightness at the orbital phase 0.4 (phase 0.0 is the 
minimum of the light curve modulation). The mass ratio 
of q = Mms/Mwd = 2 is assumed because a high mass 
transfer rate of > 1O _6 M yr" 1 is suggested (e.g. South- 
well et al. 1996). The non-irradiated surface temperature 
of the lobe-filling MS companion is assumed to be 13,000 
K (e.g., Fagotto et al. 1994). 

4. We assume the original mass transfer rate of 2 x 
10~ 6 M Q yr -1 as a reasonable value of thermally unsta- 



ble mass transfer for q = 2 and Mms = 2.6 M . This 
value is realized whenever the wind does not suppress the 
mass transfer. So, we start a rapid mass accretion to the 
WD at the rate of 2 x 10~ 6 M Q yr -1 (denoted by D-mark in 
Fig. 4). The WD envelope gradually expands and even- 
tually blows a wind (see the mid and bottom panels of 
Fig. 4). This causes a quick rise of V"-mag. Here we 
assume the disk expansion of (a, (3) = (3.0,0.05) in the 
wind phase. It takes a few to several dynamical timescales. 
We adopt 3 days. The massive winds having the rate of 
Mwi nd - 1 x 10~ 8 - 1 x 10~ 6 M Q yr" 1 certainly absorb 
supersoft X-rays (e.g., Southwell et al. 1996). 

5. This rapid mass accretion is suppressed by the strong 
wind, because the wind hits and then strips-off the very 
surface layer of the MS companion (see Hachisu et al. 
1999a, for a description of the stripping effect by winds). 
We stop the rapid mass accretion when the wind mass loss 
rate increases and the mass-stripping rate overcomes the 
mass transfer rate of 2 x 1O _6 M yr -1 . This happens 
when the WD photosphere expands to the value shown 
in Figure 4 (A-mark). However, we still keep a low rate 
(~ 1 x 1O~ 7 M yr -1 ) during the period of draining of the 
disk (a few times viscous timescale of the disk ^40 — 60 
days, e.g., Reinsch et al. 2000) even after the rapid mass 
accretion stops. 

6. When the optically thick winds stops (O-niark), we 
expect copious supersoft X-rays as observed in the U 
Sco 1999 outburst (Kahabka et al. 1999). Then the 
photosphere of the WD shrinks from 0.087 R Q (T p h = 
280,000 K) to 0.038 R Q (430,000 K) in 20 days. The 
photosphcric temperature is consistent with the blackbody 
fitted temperatures (30 — 40 eV) of the supersoft X-rays 
(e.g., Schaeidt et al. 1993) as shown in Figure 4. The 
disk shrinks to a normal size of (a, (3) = (0.8, 0.05) in a 
week and then the edge of the disk begins to flare up, i.e., 
(a, (3) = (0.8,0.1) in another week, as described in §3. 

7. Rapid mass accretion to the WD resumes (D-mark) 
~ 20 days after the wind stops. We expect that the mass 
transfer rate gradually increases just after the wind stops 
but it takes 10 — 20 days to reach the WD surface. This 
timescale may be regarded as a viscous timescale of the 
accretion disk (Reinsch et al. 2000). Then, the photo- 
sphcric radius of the white dwarf expands in 20 days and 
eventually blows winds. Then, go to 4 and repeat the 
cycle. 

Here we assume the apparent distance modulus of 
(m - M) v = 18.25(intrinsic)+0.4(absorption)= 18.65 
to RX J0513. The calculated color of (B - V) c = 
—0.26 is roughly consistent with the observational color 
of (B - V) ~ -0.14 ± 0.04 (Pakull et al. 1993), 

(B - V) 0.112 ± 0.012 (Cowley et al. 1993), or 

(B - V) 0.04 (Cowley et al. 2002) for the color 

excess of E(B - V) ~ 0.13 (Gansicke et al. 1998). 

5. DISCUSSION 

Thus, we are able to reproduce the transition between 
the optical high and low states of RX J0513 as shown 
in Figure 4. Our new model is essentially different from 
the previous expansion/contraction models of the white 
dwarf photosphere in the following sense: The emer- 
gence/decay timescale of the observed supersoft X-rays 
is as short as one or two days (Reinsch et al. 2000). 
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This short timescale cannot be explained by the expan- 
sion/contraction of the envelope. We have calculated the 
timescales of photosphcric contraction based on our WD 
envelope model (Kato & Hachisu 1994) and show them in 
Figure 5 against the WD mass. The contraction timescale 
depends on the envelope mass of the WD and it is short- 
est just after the optically thick wind stops as shown in 
Figure 4. In this phase, the timescale of contraction by a 
factor of four is about three times the Kelvin-Hclmholtz 
timescale. It is clear that the timescale of contraction by 
a factor of four cannot be shorter than six days even if the 
WD mass is very close to the Chandrasckhar mass. We 
may conclude that only the expansion/contraction model 
of the white dwarf photosphere cannot reproduce a very 
short timescale of X-ray on/off. We need another mecha- 
nism. 

On the contrary, optically thick winds completely ob- 
scure supersoft X-rays within one day after it starts. 
The wind mass loss rate increases from zero to ~ 2 x 
10~ 8 M© yr _1 in a day when it starts, or it decreases from 
~ 1 x 10~ 8 M© yr -1 to zero in a day when it ceases, as 
easily seen from Figure 4. This wind mass loss rate of 
1 x 10 _8 M© yr _1 is large enough to absorb supersoft X- 



rays (see, e.g., discussion of Southwell et al. 1996). 

In the scenario proposed by Hachisu et al. (1999a, b), 
all the progenitor systems of Type la supernovae once ex- 
perience the optically thick wind phase during their way 
to Type la supernova. RX J0513 may be an example for 
such a wind phase. We have calculated the efficiency of 
mass accretion to be about 35% for the case of Figure 4. 
The other 65% is blown in the wind or stripped off from 
the companion surface and eventually lost from the binary 
system. Assuming this efficiency being kept, we can expect 
that the white dwarf in RX J0513.9— 6951 will explode as 
a Type la supernova after about 0.6M© is transferred even 
when the WD mass is 1.2 M©. 

We thank K. Matsumoto for providing us machine read- 
able data sets of the CI Aql 2000 outburst. We are also 
greatly indebted to F. Meyer and Emi Meycr-Hofmcistcr 
for stimulating us to work on RX J0513.9— 6951 during our 
stay at the Max-Planck Institute for Astrophysics. This 
research has been supported in part by the Grant-in- Aid 
for Scientific Research (11640226) of the Japan Society for 
the Promotion of Science. 
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(a) optically thick wind phase 




(b) just after the wind stops 




(c) mass accretion makes a spray 



Fig. 1. — Configurations of our CI Aql model arc illustrated: (a) in the massive wind phase (X-ray off), (b) just after the wind stops (X-ray 
on), and (c) during a rapid mass accretion phase soon after the wind stops (X-ray on). The cool component (right) is a slightly evolved MS 
companion (1.5Mq) filling its inner critical Roche lobe. The north and south polar areas of the cool component are irradiated by the hot 
component (1.2 Mq white dwarf, left). The separation is a = 4.25Rq; the effective radii of the inner critical Roche lobes are Rjf = 1.53-Rq, 
and R2 = R2 = 1.69f?0, for the primary white dwarf and the secondary main-sequence companion, respectively, (a) In our model, we regard 
the outer edge of the disk as the transition place from optically thick to thin. When the disk surface flow collides with the companion star, 
it may run around the companion. We approximate this situation simply by the disk configuration as shown in the figure. This round flow 
hardly contributes to light curves or line profiles because its irradiation effect is much smaller than the other parts the disk, (b) The disk 
shrinks to a normal size (0.7 — 0.8 times the Roche lobe size) in several orbital periods after the wind stops, (c) A rapid mass accretion 
resumes and makes a spray around the disk edge. 
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Fig. 2. — Calculated R c light curves are plotted against the binary phase (binary phase is repeated twice from —1.0 to 1.0) together with 
the observational points (Matsumoto ct al. 2001) of the two epochs (open squares correspond to the brightness on 2000 July 17 and open 
circles are the brightness on 2000 October 1). Solid lines denote our calculated models with the inclination angle of i = 69° but dashed lines 
are of i = 74°. Various system parameters arc shown in the figure together with the best fitted parameters. 
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Fig. 3. — A calculated _R c -light curve are plotted against time (HJD 2,450,000+) together with the observational points of the CI Aql 2000 
outburst (taken from Matsumoto ct al. 2001, 2002). Small open circles indicate observational i? c -magnitudcs. The light curve connects the 
brightness at binary phase 0.4 (solid), which is roughly corresponding to the upper bound of the light curve, or at mid eclipse (binary phase 
0.0) (dashed), which gives us the lower bound. The apparent distance modulus of (m — M)r = 13.28 is used. The other system parameters 
are shown in the figure together with the best fitted parameters. Here we adopt the hydrogen content of X = 0.4 for refinement instead of 
X = 0.35 by Hachisu et al. (2002). 
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Fig. 4. — Top panel: V-magnitudc light curves of RX J0513.9— 6951 are plotted against time (JD 2,448,000+), together with the calculated 
color (B — V) c (lower portion). Small open circles are observational data with an error bar (Alcock et al. 1996). Thick solid lines correspond 
to our calculated light curve, connecting the brightness at orbital phase 0.4. Thin solid lines denote the calculated color (B — V) c - Starting 
a rapid mass accretion of 2 X 1O~ 6 M0 yr _1 to the WD at epoch of open square, we continue it until epoch of open triangle. Optically thick 
winds stop at epoch of large open circle. The maximum expansion of the WD is reached at epoch of open triangle, the radius of which is 
shown above the open triangle. Mid panel: The accretion rate to the WD (solid line) and the wind mass loss rate from the WD (dashed line) 
are plotted against time. The other symbols have the same meaning as those in the top panel. Third panel: The white dwarf radius (solid 
line) and photospheric temperature (dashed line) are plotted against time. The other symbols are the same as those in the top panel. 
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Fig. 5. — Timescale of contraction/expansion of white dwarf envelope with Z = 0.004 and X = 0.7 is plotted against the white dwarf mass. 
Here, we adopt time during which the white dwarf photosphere shrinks from 0.08 Rq to 0.02 Rq by a factor of four, which is about three 
times the Kclvin-Hclmholtz timescale of the white dwarf envelope. 



